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Abstract 

We study the equipotential surfaces around of a two particle system in 
3-d under a pairwise good potential as the one of Van der Waals. The level 
sets are completely determined by the solutions of polynomials of at most 
fourth degree that can be solved by standard algebraic methods. The 
distribution of real positive roots determines the level sets and provides a 
complete description of the map for the equipotential zones. Our methods 
can be generalized to a family of polynomials with degree multiple of 2, 
3, or 4. Numerical simulations of 2-d and 3-d pictures depicting the true 
orbits and equipotential zones are provided. 

Keywords: Pairwise good potential; solubility by radicals; Levels set; 
Equipotential Surfaces. 

1 Introduction 

The study of clusters of massless non-interactive particles under a potential 
fields has been approached, roughly speaking, under two view points: the one 
that considers the dynamics given by equations of Newton associated to the 
potential, and the one that considers static particles. The latter is related 
among others, to theoretical and experimental research in crystal formation, 
and plasticity, whereas the former is connected with research in control theory. 
Another line of research involving pairwise good potentials and control of chaos 
concerns to the dynamical analysis and control of micro cantilevers as single 
mode approximation for tackling interactions under Van der Waals potential, 
see for instance pQ. 
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In the study of the N-body problem, several pairwise good potential have 
been utilized for describing relative equilibria and central configurations, see for 
instance [3J, the methods there include classical dynamical systems as well as 
non- linear geometric control techniques [2]. 

We are interested in clusters of non-interactive particles under pairwise good 
potentials, our research points to the optimal control problems of optimal path 
planning, collision free navigation and crystal formation on equipotential zones 
of clusters. 

In this paper, we restrict ourselves to the study of the equipotential zones 
for a system of two static non interactive particles. 

In the literature, numerical methods of mathematical software is commonly 
used to generate graphics of levels set or orbits. However, there are not algebraic 
methods and techniques for determining the shape of the surfaces or curves or 
dots of the equipotential energy even for small clusters. The novelty of our 
approach is to solve the implicit function problem of the determination of the 
level sets using the algebraic methods for solving polynomials by radicals. With 
our methods, the exact solution of an orbit corresponds to roots of a polynomial 
with degree n < 4 as formulae involving its coefficients, arithmetic operations, 
and radicals. 

The paper is organized as follows, in section [2J we depict the notation of the 
problem, the general description of our approach, and our main result. Section[3J 
contains the description of the algebraic procedures. In section 0] we present a 
comparison of numerical and algebraic results. At the end in section[5l we derive 
some conclusions. 

2 Notation and Problems 

Let be R + = [0, oo). For a given d > 0, a smooth function P : R + — > R is 
well pairwise good potential if satisfies the following conditions: 

1. Infinite rejection for avoiding, destroying or collapsing particles, lim^o P(d) = 
oo. 

2. A negative basin around a minimum distance d* = argmin dg ( .cx>) P{d). 

3. Asymptotic attraction, lurid^foo P(d) = _ . 

An elementary model for the Van der Waals that complies the previous prop- 
erties is the following: B (d) = 4r — Jz- This function has its minimum at d* = 1, 
B {d*) = —1. For a group of three particles pi — (xi,yi, z,-) € R 3 , i = 1, 2, 3, with 
the Euclidean metric v/dy = d(pi,pj). We consider the problem of determin- 
ing the orbits of a system of two given particles and one free under potential, 
hereafter shall be called B 2 , that is, B 2 (px,P2,P3) = Ekkjo ^ = 
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Without lost of generality, we can assume pi = (— 1,0, 0),P2 — (5, 0,0) 
where I > is the distance for a given system of two fixed particles. Let be 
p = (x, y, z) £ K 3 a free particle, then the complete potential of this system is 
B 2 {x,y,z) = ^ (x+i _ )2+y2+z2 y - (x+J _y2 +y 2 +z 2 + ( (x _i_ )2+y2+z2 y - ( x -^ +y 2 +z 2 + 
K h 

where Ki — -k — is the corresponding potential between p\ and p 2 , which is 
constant. 

An orbit of B 2 of value G is the set {(x, y,z) £ R 3 \ B 2 (x, y, z) = G}, 
G £ [mi, 00) where I > 0, and mi = vni~n-( x ,y,z)£R a B 2 (x, y, z). 

The symmetry of B 2 (x, y, z) with respect to second and third axis allows the 
reduction of the problem to R 2 . Moreover, B 2 (x,y) — B 2 (x, — y) = B 2 {—x,—y) = 
B 2 (~x,y) . 

Therefore, it is only necessary to consider the orbits of B 2 of value G as 
the sets 0(l,G) = {(x,y) £ E+ \B 2 (x,y) = G}, I > 0, G £ [mi, 00). The 
corresponding 3-d model can be constructed by appropriate rotations. 

Our methodology is based on: (1) The algebraic methods of Cardano and 
Ferrari; (2) For a given G — G x , the equation B 2 (x,y) — G yields a third or 
fourth degree polynomial with coefficients in the ring R[x], that we solve to 
obtain an exact root r(x,G,l); (3) If r(x,G,l) > is a root, then 0(l,G) — 

I (a;, y) \ x > 0, y = \fr{x, G, T)\ . We state now the main proposition 

Proposition 1. Given a system of two particles, p\ = (— i,0, 0),£>2 = (|i0, 0). 
The orbits 0(l,G) of B 2 , for G £ [mi, 00), correspond to the positive roots of 
the third and fourth degree polynomial obtained from the equation: 

B 2 (x,y)-G = Q. (2.1) 



Proof. The equation 12.11 gives: 
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- 1? + y 2 ) + if + v 2 ) 2 - i K i + g) ((* + 1) 2 + y 2 ) 2 - 1) 2 + v 2 ) 2 )- 

It has power on x or y as 8, 6, 4, 2 when {Ki + G) 7^ and it has power on 
x or y as 6, 4, 2 when (Ki + G) = 0. 

We let now u ~ y 2 the following third and fourth degree equations are 
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obtained: 



= (l024Ga; 2 + 256G + 1280 + 1024a; 2 ) u 3 + 

(352 + 256Gx 2 + 3328a; 2 + 1536Ga; 4 + 1536a; 4 + 96G) u 2 + 

(-256Ga; 4 + 1024Gx 6 - 576a; 2 - 48 - 64Gx 2 + 2816a; 4 + 16G + 1024a; 6 ) u - 15 + 

G - 16Ga; 2 + 96Gx 4 - 256Gx 6 + 256Gx 8 - 848a; 2 - 672a; 4 + 768a; 6 + 256a; 8 , (2.2) 

= (256G + 2bQKi) u 4 + 

(l024Gx 2 + 256G + 1280 + 1024a; 2 ) u 3 + 

(352 + 256Ga; 2 + 3328x 2 + 1536Ga; 4 + 1536a; 4 + 96G) u 2 + 

(-256Ga; 4 + 1024Gx 6 - 576a; 2 - 48 - 64Ga; 2 + 2816a; 4 + 16G + 1024a; 6 ) u - 15 + 

G - 16Ga; 2 + 96Ga; 4 - 256Ga; 6 + 256Ga; 8 - 848a; 2 - 672x 4 + 768a; 6 + 256a; 8 . (2.3) 

The roots of 12.21 and 12.31 are constructed by the methods of Cardano 
and Ferrari in the complex plane C. Therefore, 0(1, G) = {(x, y) \ x > 0, y = 
y/r(x,G,l),r(x,G,l) > 0}. □ 

Remark 2. The construction of the roots is not done by mathematical software, 
but it is done by finding and replacing the parameters of the formulas of the 
methods of Cardano and Ferrari with the coefficients of the polynomials of the 
equations 12.21 and 12.31 See Appendix 1. 

Proposition 3. For the system of the previous proposition, when 1=1, we have 

1. mi = —3. 

2. 0(1, —3) = {(0, -^r)}- It is a point the root of fourth degree equation for 
Gi = -3. 

3. There are two orbits from the roots of fourth degree equation with Gi e 
[m,-l). 

4. There is one orbit from the roots of fourth degree equation with G\ G 
(-l,oo). 

Proof. It follows by direct substitution of I = 1 in the resulting roots of the 
previous proposition. □ 

Remark 4. Also for I = 1 the figure [1] depicts some examples of the orbits. 
The free point and the other two particles form an equilateral triangle of size 
1, and minimum is -62(0, ^) = —3 = mi. There is one orbit from the roots 
of third degree equation with G\ = —1, the details are given in section |H It 
is easy to prove that polynomials Ax 2k + Bx k + C, Ax 3k + Bx 2k + Cx k + D, 
Ax 4k + Bx 3k + Cx 2k + Dx k + E k > 1 are solvable by radicals. 
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3 Solving Polynomial by radicals 



The Galois theory is the algebraic framework for the study of roots of poly- 
nomials. It focuses on the construction of formulae using radicals instead of 
in numerical estimations. For numerical estimation there is for instance the 
well know method of Newton-Raphson. It is known that there is no a general 
method or formula for finding the roots of a polynomial with degree n > 5. 
Here, we are interested in applying the known methods of Cardano and Ferrari 
for polynomials of degre 3 and 4 which are solvable by radicals. Appendix 1 
depicts the formulas of Cardano and Ferrari. 

The cubic polynomial f(x) = x 3 +ax 2 +bx+c, was solved at the 16th century 
by more complicated formula found simultaneously by Ferro and Tartaglia. For 
the fourth degree polynomial f(x) — x A + ax 3 + 6a; 2 + cx + d, Ferrari provided 
a procedure. These methods were published by Cardano in the Ars Magna at 
1545. In 18th century, Lagrange unified these methods for a polynomial with 
degree n < 4 using what now is know as the resolvent of Lagrange. This method 
involves an auxiliary equation with a polynomial with degree less than one, 
i.e., it uses, for example, a third degree polynomial to solve the fourth degree 
polynomial. In fact, this is the procedure for solving the quadratic polynomial. 
However, this method fails for the polynomials of degree 5, because the auxiliary 
polynomial is of degree 6. Ruffini at 1799, and Abel at 1824 proved that there is 
not a general formula using radicals for finding the roots of quintic polynomial. 
Galois by 1832 showed how to associate to each polynomial f(x) a subgroup 
Gal(f) of the symmetric group, call the Galois group of f(x), and established 
the following result, for details see for instance [I] 

Theorem 5. (Galois) A polynomial / is soluble by radicals if and only if its 
group Gal(f) is soluble. 

As an application of these algebraic techniques of Cardano and Ferrari, we 
have the following result. 

Proposition 6. Given a system of two particles, p% = (— ^,0, 0),P2 — (^,0,0), 
the critical points of B2 (x, y) are the following 

1. (0,0) V/ > 0. 

2. (±x*, 0), x* is a positive root of the polynomial obtained from substituting 
y = in -g-Bz (x,y) with I > 2. The points (±x*,0) are collinear with 

Pl,P2- 

3. (0, ±y*) where y* is a positive root of the polynomial obtained from sub- 
stituting x = in -§zB2 (x, y). The points (0,±y*) correspond to the 
oppositive vertex of an isosceles triangle of the side with vertices p\,p2- 

Proof. The polynomial to solve comes from the first optimality condition: -§^B2 (x, y) 
= and 4-B2 (x, y) — 0, where 
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—B 2 (x,y) = -4 o H ? 

^ ( (x± / )2 + y 2)3 (^±1)2 + ^ 

-B 2 (x,y) - 



9y ((x±!) 2 + y 2 ) 3 ((^ ± + y 2 ) 2 

Note that from the previous equations, we have that VB 2 (0, 0) = 0, for all I > 0. 
Using that y = gives -§^B 2 [x, 0) = 0, the equation for -§^B 2 (u 5 , 0) = (x = u 2 
is changed) is 

-64u 3 +(64-16; 2 )w 2 + (20Z 4 + 1602 2 )u+(20Z 4 -32 6 ) = 0. (3.1) 

The positive roots of the previous third degree equation give the optimal 
points. In a similar way, using that x — gives £^B 2 (0,u^) = 0, the equation 
for -J^B 2 (0,u?) = gives 

4y 2 + I 2 - 4 = (3.2) 

The real roots of the previous second degree equation give the optimal points 
for / € [0, 2). The interval of I come from the restriction of the discriminant of 
the quadratic equation. □ 



4 Levels set of the potential B2 

The proposition [5] gives similar results for B 2 instead of LJ 12-6 to the 
triangular and collinear choreographes of the 3-body system discussed in [3]. 
We remark that our approach is done by basic analysis and algebraic techniques 
for the static case. 

Proposition 7. Given a system of two particles of the proposition [6l and I £ 
(2,0). If a free particle is at (0,y), where y = ±\\J2 2 — I 2 . Then the free 
particle is on a local point. Moreover, from the free particle's point of view, the 
two particles act as one virtual particle of double B potential when I 0. 

Proof. Without loss of generality, let y\ — \^2 2 — I 2 . The free particle is 
perpendicular to middle point of the side with vertices {pijPi)- The result follows 
immediately from proposition [B] and the equation l3.21 which corresponds to the 
condition \7B 2 (0,yi) = 0. Also, the free particle position goes to (0,2) when 
I — > 0, which is the double of the optimal distance (d* = 1) of the function 
B. □ 

Proposition 8. Given a system of two particles of the proposition [6j with 1 = 1. 
Then the minimal points of a free particle is at (0,y), where y = The 
particles are the vertices of an equilateral triangle. 
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Figure 1: The true orbits 
for 0(1,-3), 0(1,-2.75), 
0(1,-2.43755), 0(1,-2) 
0(1, -1) , 0(1, 0) , and 0(1, 100). 



Figure 2: 3-d model of the level 
sets of Bi(x,y,z) = G, where G G 
[3,oo). 



Proof. Without loss of generality, let y\ — \\fS, then the distance between the 
three particles is 1. The Hessian of the system at (0,yi) is 

|V 2 B 2 (0,^)| = -64(-3). 



The optimality follows from proposition [5] and |V 2 _B 2 (0, yi)| > 0. 



□ 



Proposition 9. Given a system of two particles of the proposition [5J and I G 
(2, 0). Then the orbit of Bi of value G y = B 2 (0, y) is the pointed set 0(/, G y )= 
{(0,y) g R+ 2 } where y = W2 2 - I 2 . 



Proof. By the symmetry the space is restricted to K + . The Hessian of the 
system at (0, y) is 

|V 2 B 2 (0, IV2 2 -Z 2 )| = -64Z 2 (l 2 - 4) . 

The polynomial — 64Z 2 (l 2 — 4) is strictly positive for I G (0, 2). This means 
that the positive root of the equation 13.21 corresponds to the optimality condi- 
tions VB 2 (0,y) = and |V 2 B 2 (0,y)| > 0. Therefore 0(1, G y ) is the minimal 
point (0, y). □ 



For G = — 3 the four roots of fourth degree equation 12.31 (with I = 1) are 

4- 



n (x) = \ - x 2 - |V4x 2 

T3 (x) 



X 



Mx + 1, r 2 (a;) = \\fAx 2 



W4a; 2 + Aix + 1, r 4 (x) 



Aix 



+ Aix + 1 



x 



Only for x — 0. r 2 (0) = 4 > 0, therefore the orbit of 0(1, —3) is the point 
set {(0, W 4)}. Which is the optimal point of the equilateral triangle of size = 1. 
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Here, our method gives the result without appealing to the optimal conditions 
as in proposition [8] and in proposition [9j 

The detailed analysis of the third and fourth degree equations of proposi- 
tion [1] for I = 1 was done for verifying the following: 

1. mi = —3. 

2. 0(1, —3) = {(0, 2 )}' i- e -' ^ * s a P om t the root of fourth degree equation 
for Gi = -3. 

3. There are two orbits from the roots of fourth degree equation with G\ € 

K-i) 

4. There is one orbit from the roots of third degree equation with G\ = — 1. 

5. There is one orbit from the roots of fourth degree equation with G\ G 
(-l,oo). 

Figure [T] depicts some examples of the orbits that comply the previous state- 
ments. Figure [2] depicts a 3-d model of the level sets for I = 1 and G <E [—3, 00). 

5 Conclusions 

To our knowledge this is the first totally complete description of the level 
sets of a system of two particles under a Van der Waals Potential. 

We have presented a novel approach for the study of orbits of systems of 
non interactive particles. We tackle the problem by reducing the analysis of the 
cquipotential zones to the one of finding the roots of certain polynomials. A 
constructive method provided by the methods of Cardano and Ferrari, yields 
a complete factorization of the polynomials and consequently an analytical de- 
scription of equipotential zones. 

Acknowledgement. Thanks to the Department of Basic Sciences at UAM, 
Azcapotzalco for supporting this research. 

Appendix 1. Formulae of Cardano and Ferrari for 12.11 

The third degree polynomial is obtained from the equation 12.11 This equation 
gives the third degree equation 12.21 with u = y 2 . Its polynomial corresponds 
to Ax 3 + Bx 2 + Cx + D = 0. And its root come from the following formula of 
Cardano: 

Xl = ((-I(A(( jB ) ( A)-i)3_I(( jB )(A)- 1 )((C)(A)- 1 )+(( J D)(A)- 1 )) + (i((C) 
(A)- 1 - l((B)(A)-i?f + I(A ((jB )(A)-i)3 _ i((B)(A)-i)((C)(A)-i) + ((D) 

(A)- 1 )) 2 )^)^)+((-I(^((B)(A)- 1 ) 3 -I((S)(A)- 1 )((C)(A)- 1 ) + (( J D)(A)- 1 ))- 
(i((C)(A)- 1 -i((B)(A)- 1 ) 2 ) 3 + I(A(( jB )(A)- 1 )3-I(( J B)(A)- 1 )((C)(A)- 1 ) + 
((D)(A)-i))*)h)i)-l((B)(A)-i) 

The fourth degree polynomial is obtained from the equation 12.11 For the 
equation 12.31 with polynomial type Ax 4 + Bx 3 + Cx 2 + Dx + E = the four 
corresponding roots come from the following four formulae of Ferrari. 
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The first root is 

(C)(A)-i) + (((-I (-^(-3(5)2(^-2 + (0(A)- 1 ) 3 + I(-f (B) 2 (A)- 2 + (C) 

(^)" 1 )(-2l6( B ) 4 (^)" 4 +^( C )( B ) 2 ^)" 3 -3( B )(^)(^)" 2 +(^)(^)" 1 )-|(l( B ) 3 
- 1{B){C){A)-* + (D)^)- 1 ) 2 ) + (i (-^(-1(5)2(^-2 + [C){A)-y + 

i(-|(S) 2 (A)- 2 + (C)(A)- 1 )(-J 5 (i?) 4 (A)- 4 + 1 3 g (C)( J B) 2 (A)- 3 -I(i?)( J D)(A)- 2 + 

(E)(A)-i) i(i(B) 3 (A)" 3 - i(B)(C)(A)- 2 + PXA)- 1 ) 2 ) 2 + ^(-^(-§(B) 2 

(A)- 2 + (C){A)^) 2 (-3is(B) 4 (A)- 4 + ^(C)(B) 2 (A)- 3 - i(i3)(D)(A)- 2 + 

(^(A)" 1 )) 3 )*))*) - l(_^(_| (B) 2 (A) -2 + (O(A)- 1 ) 2 - (-^(5)4(^-4 + 

^(C)( J B) 2 (A)- 3 -i( J B)P)(A)- 2 + (£;)( J 4)- 1 ))(((-i(- T i g (-|(i?) 2 (A)- 2 + (C) 
(A)- 1 ) 3 + I(-|(i?) 2 (A)- 2 + (C){A)^)(-^{B)\A)-± + ±{C)(BY{A)-z 
\{B){D){A)~ 2 + (E)(A)~ 1 ) I(I(i?) 3 (A)- 3 - \{B){C)(A)-* + (D)(A)~ 1 ) 2 ) 

+(i(-TM(-|( B ) 2 (^)' 2 +( c )(^)" 1 ) 3 +|(-f( s ) 2 (^)" 2 +( c )(^)" 1 )(- 2 l6( B ) 4 

(A)- 4 + ^(C)(i?) 2 (A)- 3 - i(B)P)(A)- 2 + (^(A)- 1 ) - I(I(B) 3 (A)- 3 - 1(B) 

(C) (A)-i + (D)(A)~ 1 ) 2 ) 2 + ^{-^(-l(Bf{A)- 2 + (C)(A)^) 2 (-^(S) 4 
(A)- 4 + ^(<7)(B) 2 (A)- 3 - |(B)(£>)(^)-=» + (^(A)" 1 )) 3 )*))*)" 1 ))*) - (-(3 
(-|(B) 2 (A)- 2 + (CJ(A)- 1 ) + 2(-|(-|( J B) 2 (A)- 2 + (CO(A)" 1 ) + (((-|(-^g 

( _| (jB )2 (j4) -2 + (0(A)- 1 ) 3 + |(-|(B) 2 (A)- 2 + (C)(^)- 1 )(- 5 i 5 (B) 4 (A)-' 1 + 

^(C)(i3) 2 (A)- 3 -I(i3)(I))(A)- 2 + (i?)(A)- 1 )-i(I( J B) 3 (A)- 3 -I( J B)(C)(yl)- 2 + 

(D) (A)- 4 ) 2 ) + (i(- i (-|(B) 2 (A)- 2 + (C)(A)^f + I(-|(S) 2 (A)- 2 + (C) 

(^)" 1 )(-2l6( B ) 4 (^)^+il(^)( B ) 2 ^)" 3 -K B )(^)(^)" 2 +(^)(^)" 1 )-i(5( B ) 3 
(A)- 3 -i(i?)(C)(A)- 2 + ( J D)(A)- 1 ) 2 ) 2 + ^(-i(-|(i?) 2 (A)- 2 + (C)(A)- 1 ) 2 - 

(-4(B) 4 (A)- 4 + i(C)(B) 2 (A)- 3 - I(B)P)(A)- 2 + (^(A)- 1 )) 3 )!))!) - 

(D)(A)- 2 + (£;)(A)- 1 ))(((-I(- T 1 ? (-|(B) 2 (A)- 2 + (C)(A)-y + §(-|(£) 2 
(A)- 2 + (0{A)-i)(-MB)*{A)-* + 1(C)(B) 2 (A)- 3 - i(B)(£>)(A)- 2 + (25) 
(A)- 1 ) - I(I(i3) 3 (A)- 3 - i(S)(C)(A)- 2 + (^(A)- 1 ) 2 ) + (I(-J-(-§ (B ) 2 

(A) - 2 + (C)(A)- 4 ) 3 + l(-l(B) 2 (A)- 2 + (C)(A)-i)(-MBf(A)-* + ±(C) 

(B) 2 (A)-z - \{B)(D){A)- 2 + (E)(A)-i) I(I(S) 3 (A)" § - \(B){C)(A)- 2 + 

(D) (A)-i) 2 ) 2 + 2 f (-±(-l(B) 2 (A)- 2 + (C)(A)-i) 2 - (-JL(B) V)~ 4 + ^ {C) 

(B) 2 (A)- 3 -I(i?)( J D)( J 4)- 2 + (ii;)(A)- 1 )) 3 )i))l)- 1 ) + 2(I(i?) 3 (A)- 3 -I(i?)(C) 
(A)- 2 + (D)(A)-i)(((-l(B) 2 (A)- 2 + (C)(A)-i) + 2(-l(-l(B) 2 (A)- 2 + (C) 
(A)- 1 )+(((-i(- T i I (-|(i?) 2 (A)- 2 +(C)(A)- 1 ) 3 +I(-f(i?) 2 (A)- 2 +(C)(A)- 1 ) 
{-£-{B)\A)-± + ±(C)(B) 2 (A)- 3 \{B){D){A)- 2 + (E)(A)~ 1 ) I(±(£) 3 
(A)-* - ±{B){C)(A)- 2 + {D)(A)-^) 2 ) + ( i (- i (-|(S) 2 (A)- 2 + (C)(A)-if + 
i(-|(B) 2 ( J 4)- 2 + (C)(A)- 1 )(-^(i?) 4 (A)- 4 +l(C)(B) 2 (A)- 3 -I( J B)( J D)(A)- 2 + 

(E) (A)-i) i(i(B) 3 (A)" 3 - i(B)(C)(A)- 2 + PXA)- 1 ) 2 ) 2 + ^(-^(-f^) 2 
(A)- 2 + (C)(A)- 4 ) 2 - (-^(B)4(A)-4 + i(C)(B) 2 (A)" 3 - \(B){D)(A)- 2 + 
(E)(A)-^))i) l(-±(-l(B) 2 (A)- 2 + (C)(A)-i) 2 (-^(5)4(^-4 + 
i(C)( J B) 2 (A)- 3 -i(B)P)(A)- 2 + (£;)(A)- 1 ))(((-i(- T i g (-|(i?) 2 (A)- 2 + (C) 
(A)" 1 ) 3 + !(_§(B) 2 (A)- 2 + (<7)(A)-i)( 3 (B^H + ±(C)(B) 2 (A)- 3 
\{B){D){A)- 2 + (E)(A)^ 1 ) - l(l{Bf{A)^ - i(B)(C)(A)- 2 + (D)(A)^ 1 ) 2 ) 
+ (l(-T58(-|(5) 2 (A)- 2 +(C)(A)- 1 ) 3 + i(-|( J B) 2 (A)- 2 + (C)(A)- 1 )(- 5 fg( J B) 4 
(A)" 4 + l(C)(i?) 2 (A)- 3 - i(S)P)(A)- 2 + (^(A)- 1 ) - I(I(B) 3 (A)- 3 - 1(B) 

(C) (A)- 2 + W(A)- 1 ) 2 ) 2 + ±(-±(-l(B) 2 (A)- 2 + (C)(A)-i) 2 (-^(S) 4 
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(A)- 4 + ±(C)(B)*(A)-* - \{B){D)(A)-* + {E){A)-^)f)h))\)-i))h)-i))h). 
The second root is 

x 2 = -^B)(A)^ + ^(((-l(B)^A)- 2 + (C)(A)^) + 2(-l^l(Bf(A)- 2 + 
(C)(A)-i) + (((-|(- JL(-f(B)2(A)- 2 + (C)^)- 1 ) 3 + I(-|(i?)2(A)- 2 + (C) 

(^)" 1 )(-2i6(5) 4 (^)- 4 +^(c')(^) 2 (^)~ 3 -i(5)(^)(^)" 2 +(^)(^)" 1 )-i(i(5) 3 

(A)- 3 - i(5)(C)(A)- 2 + (D)^)" 1 ) 2 ) + (I (-^(-KS) 2 ^)" 2 + (C)(A) -1 ) 3 + 
i(-|(S) 2 (A)- 2 + (C)(yl)- 1 )(-J g (i?) 4 (A)- 4 +i(C)( J B) 2 (A)- 3 -I( J B)( J D)(A)- 2 + 
(^(A)" 1 ) - i(i(S) 3 (^)- 3 - i(5)(C)(A)- 2 + (^(A)- 1 ) 2 ) 2 + ^(-^(-§ (j B) 2 
(A)- 2 + (C)(A) -1 ) 2 - (-^(B) 4 (A)- 4 + jq(C)(B) 2 (A)~ 3 - \{B){D){A)-* + 
{E){A)^)f )k))\) i(-^(-3 (B) 2 ( ^-2 + (^(A)- 1 ) 2 - (-^(B) 4 ^)- 4 + 
^(C)( J B) 2 (A)- 3 -i(B)P)(^)- 2 + (ii;)(A)- 1 ))(((-i(- T i g (-|(i?) 2 (A)- 2 + (C) 
(A)- 1 ) 3 + I(-|(i?)2(A)- 2 + (C)^)- 1 )^^^) 4 ^)- 4 + ^(C)^) 2 ^)- 3 - 
\{B){D)(A)-* + (E)(A)~ 1 ) - |(I(i3) 3 (A)- 3 - i(B)(C)(A)- 2 + (B)^)- 1 ) 2 ) 

+ (3(-TM(-|( B ) 2 (^)" 2 +(^)(^)" 1 ) 3 + l(-f( 5 ) 2 (^)" 2 + (G)(^)" 1 )(-2f6( B ) 4 
(A)- 4 + l(C)(i?) 2 (A)- 3 - i(B)(D)(A)- 2 + (E)(A)-i) I(I(B) 3 (^)- 3 - 1(B) 

(C) (A)- 2 + (^(A)- 1 ) 2 ) 2 + ^(-^(-KB)^^)- 2 + (C)(A)^f - (-^(S) 4 
(A)- 4 + i(C)(i?) 2 (A)- 3 - i(i?)(D)(A)- 2 + (B)^)- 1 )) 3 )^))*)- 1 ))^) + (-(3 
(-f(i?) 2 (A)- 2 + (CJ(A)- 1 ) + 2(-f(-f( J B) 2 (A)- 2 + (CO(A)" 1 ) + ({{-\(-^ 
(-|(i?) 2 (A)- 2 + (C)^)" 1 ) 3 + I(-§(B)2 (j4 )- 2 + (^(^-^(-^(i?) 4 ^)- 4 + 
^(C)(i3) 2 (A)- 3 -I(i3)(D)(A)- 2 + ( J B)(A)- 1 )-i(I( J B) 3 (A)- 3 -I( J B)(C)(yl)- 2 + 

(D) (A)-^) + (l(-^l( B nA)- 2 + (C)(A)-if + |(-|(S) 2 (A)- 2 + (C) 

(^)" 1 )(-^( B ) 4 (^)^+^(^)( B ) 2 ^)" 3 -3( B )(^)(^)" 2 +(^)(^)" 1 )-|(|( B ) 3 
(A)- 3 -i(i?)(C)(A)- 2 + ( J D)(A)- 1 ) 2 ) 2 + i(- 1 L(-|(i?)2(A)- 2 + (C)(A)- 1 ) 2 - 

(-^(B) 4 (A)- 4 + ^(^(B) 2 ^)- 3 - \{B){D){A)-* + (^(A)- 1 )) 3 )^))!) - 
^-^(-IW'^^ + ^^^^'-l-Je^^^ + ll^)^) 2 ^)" 3 -!^) 
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(D)( 1 l)- 2 "+ (E)(A)-i))(((-U-^(-({B) 2 (A)-2 + (C)(A)-y + i(-|(B) 2 
(A)- 2 + (CX^X-^B)^)-* + ^(C)(B) 2 (A)- 3 - i(B)(£>)(^l)- 2 + (£) 

(A) - 2 + (C)^)" 1 ) 3 + I(-|(S) 2 (A)- 2 + (C)(A) _1 )(— (_B) 4 ( J 4) -4 + i(C) 

(B) 2 (A)- 3 - i(B)P)(A)" 2 + (E)(A)~ 1 ) - I(I(B) 3 (A)-§ - \{B){C){A)-* + 
(^(A)- 1 ) 2 ) 2 + ^(-^(-K^^)- 2 + (C)(A)^f - (-^(B) 4 (A)- 4 + ±(C) 
( J B) 2 (A)- 3 -I(i?)( J D)( J 4)- 2 + (ii;)(A)- 1 )) 3 )^))^)- 1 ) + 2(I(i?) 3 (A)- 3 -I(i?)(C) 
(A)- 2 + (^(^-^(((-KB) 2 ^)- 2 + (C)(A)-') + 2(-|(-|(B) 2 ( J 4)- 2 + (C) 
(A)- 1 ) + (((-i(- T ^(-|(i?) 2 (A)- 2 + (C)(A)- 1 ) 3 + I(-|( J B) 2 (A)- 2 + (C)(A)- 1 ) 
(-4(i?) 4 (A)" 4 + ^(C)(B) 2 (A)- 3 - i(B)P)(A)- 2 + (^(A)- 1 ) - I(i(B) 3 

(-^) - f (sxcx^r 2 + WW- 1 ) 2 ) + (i(-TM (-|(s) 2 (^)- 2 + (c)^)- 1 ) 3 + 

i(-|(B) 2 (A)- 2 + (C)(yl)- 1 )(-^(i?) 4 (A)- 4 +l(C)( J B) 2 (A)- 3 -I( J B)( J D)(A)- 2 + 
(^(A)" 1 ) - I(I(S) 3 (yl)- 3 - I(B)(C)(A)- 2 + (D)^)" 1 ) 2 ) 2 + ±(-±(-l(Bf 
(A)- 2 + (C)(A) -1 ) 2 - (-^(B) 4 (A)- 4 + i(C)(B) 2 (A)" 3 - \(B){D)(A)-* + 

{E){A)^)f)^) |(-^(-|(B) 2 (^)- 2 + (CXA)- 1 ) 2 - (-^(B) 4 ^)- 4 + 
i(C)( J B) 2 (A)- 3 -i(B)P)(A)- 2 + (£;)( J 4)- 1 ))(((-i(- T i I (-|(i?) 2 (A)- 2 + (C) 
(A)- 1 ) 3 + I(-|(B)2(A)- 2 + (C)( J 4) _1 )(— 2fg(B) 4 (A) -4 + i(C)(B) 2 (A)- 3 - 
i(B)p)(A)- 2 + (^(A)- 1 ) - I(I(i?) 3 (A)- 3 - i(B)(C)(A)- 2 + (D)(A)-i f) 

+(i(-TM(-|( B ) 2 (^)" 2 +( c )(^)^ 1 ) 3 +i(-f( B ) 2 (^)" 2 +( c )(^)" 1 )(-2i6( B ) 4 

(A)- 4 + 1(CXB) 2 (A)- 3 - |(S)P)(A)- 2 + (^(A)- 1 ) - I(I(B) 3 (^)- 3 -1(B) 
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(A)- 4 + ±.(C)(BnA)-z \{B){D){A)-z + (^(A)- 1 )^)^)^)-!))^)-!))^). 
The third root is 

x 3 = -I( J B)( J 4)- 1 +I(-(((-|( J B) 2 (A)- 2 + (C)(A)- 1 )+2(-|(-|( J B) 2 ( J 4)- 2 + 
(C)(A)-') + (((-^(-^(-KB) 2 ^)- 2 + (C)(A)-^ + I(-|(i?) 2 (A)- 2 + (C) 
(^)" 1 )(-j6(5) 4 (^)- 4 +^(C')( J B) 2 (A)- 3 -i(i?)( J D)(A)- 2 +(£;)(A)- 1 )-i(i(i3) 3 
(A)- 3 - i(B)(C)(A)- 2 + (D)(A)-if) + ( I (-^(-KS) 2 ^)" 2 + (CXA)- 1 ) 3 + 
i(-|(S) 2 (A)- 2 + (C)(A)- 1 )(-J 5 (i?) 4 (A)- 4 + 1 3 H (C)( J B) 2 (A)- 3 -I( J B)( J D)(A)- 2 + 
(^(A)" 1 ) - I(I(B) 3 (yl)- 3 - I(B)(C)(A)- 2 + PXA)- 1 ) 2 ) 2 + ±(-±(-l(B f 
(A)- 2 + (C){A)^f (--^ {B )\A)-± + ±(C)(B)*(A)-z ~ \{B){D){A)-* + 
{E){A)^)f)k))\) i(-^(-| (j B)2 (j 4)- 2 + (C)(A)-^ (-^(S) 4 ^)- 4 + 
^(C)( J B) 2 (A)- 3 -i(i3)P)(A)- 2 + (£;)(A)- 1 ))(((-i(- T ^(-|(i?) 2 (A)- 2 + (C) 
(A)- 1 ) 3 + I(-|(i?) 2 (A)- 2 + (C)^)- 1 )^^^) 4 ^)- 4 + ^(CXB) 2 ^)- 3 - 
\{B){D)(A)-* + (E)(A)~ 1 ) - I(I(i?) 3 (A)- 3 - \{B){C)(A)-* + (D)(A)~ 1 ) 2 ) 

+ (3(-TM(-|( B ) 2 (^)" 2 +( C )(^)" 1 ) 3 + i(-|( B ) 2 (^)" 2 + ( C )(^)" 1 )(-2i6( B ) 4 
(A)^ + l(C)(i?) 2 (A)- 3 - i(B)P)(A)- 2 + (^(A)" 1 ) - I(I(S) 3 (A)" 3 - 

(C) (A)- 2 + (^(A)- 1 ) 2 ) 2 + ^(-i^KB) 2 ^)- 2 + (C)^)- 1 ) 2 - (-^(S) 4 
(A)- 4 + i(C)(i?) 2 (A)- 3 - \{B)(D){A)- 2 + (^(A)- 1 )) 3 )^))*)- 1 ))^) - (-(3 
(-f^) 2 ^)" 2 + (CJ(A)- 1 ) + 2(-|(-|(i?) 2 (A)- 2 + (CO(A)" 1 ) + (((-|(-Tfe 
(-f^) 2 ^)- 2 + (C)^)" 1 ) 3 + K-KB) 2 ^)- 2 + (C)(A)-i)(-^.(B)\A)-* + 
i(C)(B) 2 (A)- 3 -I(B)(I))(A)- 2 + (i?)(A)- 1 )-i(I( J B) 3 (A)- 3 -l( J B)(C)(A)- 2 + 

(D) ^)- 1 ) 2 ) + J-(_f(5) 2 (yL)- 2 + (O(A)- 1 ) 3 + I(-|(S) 2 (^)- 2 + (C) 
(^)" 1 )(-j6(W)- f +^(C)(i3) 2 (A)- 3 -I(i?)(Z?)(A)- 2 +(^)(A)- 1 )-I(i(i3) 3 
(A)- 3 - I(B)(C)(,4)- 2 + (^(A)- 1 ) 2 ) 2 + ^(-^(-KB) 2 ^)- 2 + (O(A)- 1 ) 2 - 
(-^(S) 4 (A)- 4 + i(C)(B) 2 (A)- 3 - 1(B)(0)(A)- 2 + (^(A)- 1 )) 3 )^))^ - 

5(-i(-|( i? ) 2 ^)" 2 +( c )(^)" 1 ) 2 -(-2i6( i? ) 4 (^)" 4 +ii( c )( B ) 2 (^)" 3 -3( B ) 

(D) (A)- 2 + (£;)(A)- 1 ))(((-I(- T 1 ? (-|(B) 2 (A)- 2 + (C)(A)-^ + l 3 (-l(B) 2 
(A)- 2 + (C)(A)-')(-^(Bf(A)- i + ^(C)(B) 2 (A)- 3 - I(B)(£>)(A)- 2 + (25) 
(A)- 1 ) - i(i(B) 3 (A)- 3 - i(5)(C)(A)- 2 + (^(A)- 1 ) 2 ) + (I(-J-(-f(B) 2 

(A) - 2 + (C)^)- 1 ) 3 + ±{-l(B)HA)- 2 + (C)(A)^)(-^-(Bf(A)^ + ±(Q 

(B) 2 (A)- 3 - I(i?)(£>)(A)- 2 + (^(A)- 1 ) - i(|(B) 3 (A)- S - i(B)(C)(A)- 2 + 
ODXA)- 1 ) 2 ) 2 + i(-i(-|( B ) 2 (A)- 2 + (C)^)- 1 ) 2 - (-^(B) 4 ^)- 4 + -L(C) 

(S) 2 (A)- 3 -i(S)(£>)(A)- 2 + ( J B)(A)- 1 )) 3 )i))l)- 1 )-2(i(S) 3 (A)- 3 -i(S)(C) 
(A)- 2 + (D)(A)-i)(((-l(BnA)-* + (C)(A)-i) + 2(-|(-|(B) 2 (A)- 2 + (C) 
(A)- 1 ) + (((-i(- T i g (-|(i?) 2 (A)- 2 + (C)(A)- 1 ) 3 + I(-|( J B) 2 (A)- 2 + (C)(A)- 1 ) 
(-4(B) 4 (A)- 4 + i(C)(i?) 2 (A)- 3 - i(B)(D)(A)- 2 + (^(A)- 1 ) - I(i(B) 3 
- I(B)(C)(A)- 2 + (^(A)" 1 ) 2 ) + (i (- i (_|(B)2(A)- 2 + (C)^)" 1 ) 3 + 
i(-|(B) 2 (A)- 2 + (C)(A)- 1 )(-J n (i?) 4 (A)- 4 +l(C)( J B) 2 (A)- 3 -I( J B)( J D)(A)- 2 + 

(E) ^)- 1 ) - i(i(S) 3 (A)" 3 - i(B)(C)(A)- 2 + (.D)(A) -1 ) 2 ) 2 + ^(-^(-f^) 2 
(A)- 2 + (C){A)^f (-^ (jB ) 4 (A)- 4 + i(C)(B) 2 (A)" 3 - \(B){D)(A)-* + 
(E)(A)-^))i) l(-±(-l(Bf(A)- 2 + (C)(A)-if (-^(5)4(^-4 + 
^(C)( J B) 2 (A)- 3 -i(B)P)(A)- 2 + (£;)( J 4)- 1 ))(((-i(- T i 5 (-|(i?) 2 (A)- 2 + (C) 
(A)- 1 ) 3 + I(-|(B) 2 (A)- 2 + (C)^)- 1 )^^^) 4 ^)- 4 + i(C)(B) 2 (A)" 3 - 
i(5)P)(A)- 2 + (^(A)- 1 ) - I(I(i?) 3 (A)~ 3 - i(B)(C)(A)- 2 + (D)(A)^ 1 ) 2 ) 

+(i(-Tk(-!( B ) 2 ^)" 2 +( c )(^)" 1 ) 3 +i(-f( s ) 2 ^)" 2 +( c )^)" 1 )(-2i6( B ) 4 
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(A)- 4 + i(C)CB) 2 (A)- 3 - i(B)OD)(A)- 2 + (E)(A)^) - l(l(Bf(A)- 3 -\{B) 
(C)(A)-i + (D){A)^ff + ^(-^(-|(B)2(^)-2 + (C)(A)-^Y (-M B ) 4 
(A)- 4 + ±(C)(BnA)-z \{B){D){A)-z + (^(A)- 1 )) 3 )^))!)- 1 ))^)- 1 ))!). 
Finally, the fourth root is 

Xi = -I( J B)( J 4)- 1 +I(-(((-|( J B) 2 (A)- 2 + (C)(A)- 1 )+2(-|(-|( J B) 2 ( J 4)- 2 + 
(C)(A)-') + (((- \ (-^(-KS) 2 ^)" 2 + (C)^)- 1 ) 3 + I(-f(B) 2 (A)- 2 + (C) 

(^)" 1 )(-^( B ) 4 (^)" 4 +ra( c )( 5 ) 2 (^)" 3 -i( B )(^)(^)" 2 +(^)(^)" 1 )-|(i( B ) 3 

(A)- 3 - i(5)(C)(A)- 2 + (^(A)- 1 ) 2 ) + (I (-^(-KS) 2 ^)" 2 + (C)(A) -1 ) 3 + 
i(-|(S) 2 (A)- 2 + (C)(yl)- 1 )(-J g (i?) 4 (A)- 4 + 1 i g (C)( J B) 2 (A)- 3 -I(i3)( J D)(A)- 2 + 
(^(A)" 1 ) - |(i(B) 3 (^)- 3 - i(5)(C)(A)- 2 + (^(A)- 1 ) 2 ) 2 + ^(-^(-f^) 2 
(A)- 2 + (CXA)- 1 ) 2 - (-^(B) 4 (A)- 4 + jq(C)(B) 2 (A)~ 3 - \{B){D){A)-* + 
(^(A)- 1 )) 3 )!))!) - L(-±(-l( B nA)-2 + (C)(A)-'f (-^(5)4(^-4 + 
i(C')( J B) 2 (A)- 3 -i( J B)P)(A)- 2 + (£;)(A)- 1 ))(((-i(- T i g (-|(i?) 2 (A)- 2 + (C) 
(A)- 1 ) 3 + I(-|(i?) 2 (A)- 2 + (C)^)- 1 )^^^) 4 ^)- 4 + ^(C)^) 2 ^)- 3 - 
i(B)P)(A)- 2 + (^(A)- 1 ) - I(I(i?) 3 (A)- 3 - i(B)(C)(A)- 2 + (C)^)- 1 ) 2 ) 

+(i(-TM(-|( B ) 2 (^)" 2 +( c )(^)" 1 ) 3 +i(-|( B ) 2 (^)^ 2 +( c )(^)" 1 )(-2i6( B ) 4 

(A)- 4 + J(C)(i?) 2 (A)- 3 - i(B)(D)(A)- 2 + (^(A)" 1 ) - I(I(B) 3 (A)- 3 -1(B) 

(C) (A)- 2 + (^(A)- 1 ) 2 ) 2 + ^(-^(-f(B) 2 (A)- 2 + (C)(A)^) 2 - (-Je(B) 4 
(A)- 4 + i(C)(i?) 2 (A)- 3 - i(B)(D)(A)- 2 + (E)(A)^)r)i))i)-^) + (-(3 
(-|(B) 2 (A)- 2 + (C)(A)^) + 2(-|(-|(i?) 2 (A)- 2 + (C)(A)-i) + (((-|(-T5s 
(-f^) 2 ^)- 2 + (CXA)- 1 ) 3 + I(-|(B) 2 (,4)- 2 + (^(^-^(-^(i?) 4 ^)- 4 + 
^(C)(B) 2 (A)- 3 -I(i?)(D)(A)- 2 + (i?)(A)- 1 )-i(I(i?) 3 (A)- 3 4( J B)(C)(A)- 2 + 

(D) ^)- 1 ) 2 ) + (-|(B)2(A)- 2 + (C)(A)-if + I(-|(B) 2 (yl)- 2 + (C) 

(A)- 3 -I(i?)(C)(A)- 2 + ( J D)(A)- 1 ) 2 ) 2 + i(-^(-|(i?) 2 (A)- 2 + (C)(A)- 1 ) 2 - 
(-^(B) 4 (A)- 4 + i(C)(B) 2 (A)" 3 - i(B)(£»)(A)- 2 + (^(A)- 1 )) 3 )^))!) - 

(D)(A)- 2 + ( £ ;)(A)- 1 ))(((-I(- T i F (-f(B) 2 (A)- 2 + (C)(A)-^ + i(-|(B) 2 
(A)- 2 + (OiAr^i-MBfiA)-* + 1(C)(B) 2 (A)" 3 - i(B)(£>)(A)" 2 + (£) 
(A)- 1 ) - I(I(B) 3 (A)- 3 - i(B)(C)(A)- 2 + (^(A)- 1 ) 2 ) + i (_|(B) 2 

(A) - 2 + (C)^)- 1 ) 3 + §(-g (B) 2 (A)- 2 + 3 + i(C) 

(B) 2 (A)- 3 - i(B)(£»)(A)" 2 + (E)(A)-i) - I(I(B) 3 (A)- § - i(B)(C)(A)- 2 + 

(D) (A)- 4 ) 2 ) 2 + i(-i(-|( B ) 2 (A)- 2 + (C)(A)^f - (- JL(B) 4 (A)- 4 + !_ {c) 
( J B) 2 (A)- 3 -I(i?)( J D)( J 4)- 2 + (ii;)(A)- 1 )) 3 )^))^)- 1 )-2(I(i?) 3 (A)- 3 -I(i?)(C) 
(A)- 2 + (^(^-^(((-KS) 2 ^)" 2 + (C)(A)-i) + 2(-|(-| (B) 2 (A)- 2 + (C) 
(A)- 1 ) + (((-i(- T i g (-|(B) 2 (A)- 2 + (C)(A)- 1 ) 3 + I(-|( J B) 2 (A)- 2 + (C)(A)- 1 ) 
(- JL(B) 4 (A)- 4 + i(C)(5) 2 (A)- 3 - i(B)P)(A)- 2 + (^(A)- 1 ) - |(i(B) 3 
(A) - i(S)(C)(A)- 2 + (D)(A)-if) + (l(-^(-l(Bf(A)- 2 + (C)(A)-if + 
i(-|(B) 2 (A)- 2 + (C)(A)- 1 )(-J n (i?) 4 (A)- 4 +i(C)( J B) 2 (A)- 3 -I( J B)(I))(A)- 2 + 

(E) (A)-i) Kl(B) 3 (A)- 3 \{B){C){A)-* + {D)(A)^ff + ^(-^(-f^) 2 
(A)- 2 + (CXA)- 1 ) 2 - (-^(B) 4 (A)- 4 + ±(C)(Bf(A)- 3 - i(B)P)(A)- 2 + 
(^(A)" 1 )) 3 )^)!) - |(-^(-|(B) 2 (^)- 2 + (C)(A)-^ (-^(B) 4 ^)- 4 + 
^(C)(i3) 2 (A)- 3 -i(B)P)(A)- 2 + (^)( J 4)- 1 ))(((-i(- T i g (-|(i?) 2 (A)- 2 + (C) 
(A)- 1 ) 3 + I(-|(B) 2 (A)- 2 + (C){A)^)(-^{B)\A)-± + ±(C)(B) 2 (A)- S 
\{B){D){A)-z + (E)(A)~ 1 ) ~ |(I(i?) 3 (A)- 3 - \{B){C){A)-z + (-D)(A) -1 ) 2 ) 
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MJ(-^h(-UB)^A)-H(C)(A)-y^(-|(B)^A)-H(C)(A)- ^ )(-^(B^ 
(A)- 4 + i(C)(i?) 2 (A)-3 _ i mD){A) -2 + {E )(A)- 1 ) - i (|(B)3(A)-3 -1(B) 
(G)(A)- 2 + (^(A)- 1 ) 2 ) 2 + ^(-^(-|(B)2 (j 4)- 2 + (C)(A)-^ - {-M B )' 
(A)~ 4 + ^(C)(B) 2 (A)"3 - I(B)P)(A)- 2 + (l^)- 1 )) 3 )*))*)- 1 ))*)- 1 ))*). 
Remark 10. The construction of the roots is done by a text editor, replacing 
the parameters (A,B,C,D, and £7), of the previous root formulae with the co- 
efficients of given polynomials. The parenthesis are necessary to avoid errors 
when complex term are replaced, by example, replacing — 256Gx 4 + 1024Ga; 6 — 
576x 2 - 48 - 64Ga; 2 + 2816a; 4 + 16G + 1024a; 6 by D for the formulae of Ferrari. 
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